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By  Robert  R.  Long 
Department  of  Earth  Sciences 
The  Johns  Hopkins  University 
Baltimore,  Maryland  21218 

ABSTRACT 

A theory  is  developed  for  turbulence  in  a stably  stratified  fluid,  for  example 
in  the  experiments  of  Rouse  and  Dodu  and  of  Turner.  In  these  there  is  no  shear 
and  the  turbulence  is  induced  by  a source  of  energy  near  the  lower  boundary  of 
the  fluid.  A growing  mixed  layer  of  thickness  D appears  in  the  lower  portion  of 
the  fluid,  separated  from  the  non -turbulent  fluid  above,  in  which  the  mean  buoyancy 
gradient  is  given,  by  an  interfacial  layer  (IL)  of  thickness  h.  The  lower  mixed 
layer  has  a very  weak  buoyancy  gradient  and  the  large  buoyancy  difference  across 
IL  is  fib.  As  indicated  by  the  experiments  of  Thompson  and  Turner  and  Hopfinger 
and  Toly,  and  derived  by  the  author  in  a recent  paper,  if  u is  the  rms  horizontal 
velocity  and  JL  is  the  integral  length  scale,  u£  = K is  a constant  in  a layer  of  homo- 
geneous fluid  agitated  by  a grid  and  at  some  distance  from  the  grid.  When  there 
is  stratification,  the  fluid  motion  is  unaffected  by  buoyancy  forces  in  the  mixed 
layer  so  that  u l should  also  be  constant  in  the  lower  portions  of  the  mixed  layer. 
Since  l is  proportional  to  distance,  we  may  conveniently  suppose  that  the  source 
of  the  disturbances  is  at  a level  z = 0 where  u is  infinite  in  accordance  with  uz  = K. 
Thus  we  may  take  K to  be  a fundamental  parameter  characterizing  the  turbulent 
energy  source.  Then  z is  distance  above  the  plane  of  the  virtual  energy  source. 

If  the  upper  fluid  is  of  uniform  buoyancy,  Dfib  = v2  may  be  shown  to  be  constant 
if  we  accept  the  experimental  observation  that  h is  proportional  to  D.  In  general 
v may  be  taken  to  be  a fundamental  parameter  expressing  the  stability.  The 
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quantity  Ri  = vsDs/Ka  is  the  most  fundamental  of  the  several  Richardson 
numbers  that  have  been  introduced  in  this  problem  because,  with  its  use, 

"constants"  of  proportionality  do  not  depend  on  the  molecular  coefficients  of 
viscosity  or  diffusion  (for  high  Reynolds  number  turbulence)  or  on  the  geometry 
of  the  grid. 

_ 7 

The  theory  contains  a number  of  results:  (1)  u,D/K  ~ Ri  r,  where  u, 

2 

is  entrainment  velocity,  i.  e.  the  rate  of  increase  of  D.  This  implies  DaP 

i 

for  a homogeneous  upper  fluid  and  D«tff  for  the  upper  linear  density  field. 

The  entrainment  law  compares  with  a -f  law  suggested  by  several  experimenters; 

- -2 

(2)The  turbulent  field  in  IL  is  intermittent  with  intermittency  factor  I ~ Ri  . 

The  turbulent  patches  have  dimension  6 ~ DRi  4 ; (3)  If  the  rms  horizontal  velocity 
at  some  fixed  level  not  far  from  the  grid  is  denoted  by  ux , we  find  that  the  rms 

* -i 

velocity  near  the  interface  is  u2  ~ ux  Ri  4 ; (4)  The  buoyancy  flux  near  the  inter- 
face q0  may  be  expressed  as  q0  ~u|/D  as  suggested  by  the  author  in  an  earlier 
paper.  j 

A recent  experimental  paper  by  Hopfinger  and  Toly  proposes  that  ux  ~us. 

The  present  theory  suggests  that  the  measurements  of  the  rms  velocity  were 
made  too  far  from  the  interface. 

1.  Introduction. 

Experimental  observations  beginning  with  those  by  Rouse  and  Dodu  (1955) 
show  that  if  a stably  stratified  fluid  is  agitated,  say  at  the  bottom  of  a container, 
a mixed  layer  develops  near  the  bottom  of  depth  D increasing  with  time.  The 
mean  buoyancy  profile  observed  is  shown  schematically  in  figure  1.  The  turbulence 
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dies  out  across  an  interfacial  layer  (IL)  of  thickness  h.  The  problem  has 
great  importance  in  meteorology  and  oceanography  as  discussed  by  many 
authors,  for  example  Kraus  and  Turner  (1967). 

Cromwell  (1960)  constructed  a similar  experiment  to  simulate  the 
pycnocline  but  recent  interest  began  with  the  careful  measurements  by  Turner 
(1968).  Subsequently  many  others  have  reported  on  identical  or  similar  ex- 
periments (Brush,  1970,  Wolanski,  1972,  Linden,  1973,  Crapper  and  Linden, 
1974,  Linden,  1975,  Wolanski  and  Brush,  1975,  Thompson  and  Turner,  1975, 
Hopfinger  and  Toly,  1976).  Other  experiments  involving  shear  have  been  run 
by  Kato  and  Phillips  (1969),  Moore  and  Long  (1971),  Wu  (1973),  Kantha,  1975, 
Kantha,  Phillips  and  Azad  (1977).  Typical  of  these  is  that  of  Kato  and  Phillips 
shown  schematically  in  figure  2. 

A number  of  suggestions  have  been  made  based  on  observation.  In  the 
experiment  without  shear,  experimenters  agree  that  the  entrainment  velocity 
u,  is  given  by1 


where  Ri*  is  the  overall  Richardson  number  defined  in  the  legend  of  figure  1. 
In  addition  there  is  strong  evidence,  presented  by  Long  (1973),  Crapper  and 
Linden  (1974),  and  more  recently  by  Hopfinger  and  Toly  (1976),  that  h = aD 
where  a is  independent  of  the  Richardson  number.  In  Crapper  and  Linden's 
measurements  h was  constant  as  the  Richardson  number  varied  from  4 to  6000 


JIn  this  paper,  we  assume  that  the  overall  Richardson  number  Ri*  or 
similar  non-dimensional  number  is-  large.  If  two  non-dimensional  numbers 
A,  B have  a ratio  A/B  tending  to  a finite  nonzero  constant  as  Ri*  - 00 , we  say 
that  A is  of  order  B and  write  it  A ~B.  We  use  the  proportionality  symbol 
connecting  two  dimensional  quantities  that  vary  together. 
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so  there  can  be  little  doubt  that  a is  independent  of  Ri.  It  is  possible,  however, 
that  a may  vary  with  the  density  gradient  in  the  upper  layer. 

In  shearing  experiments  there  is  some  indication  that  u,  /u+  ~ Ri*-1  but 
considerable  doubt  has  been  raised  by  recent  experiments  (Kantha,  1975,  Kantha, 
et  al,  1977)  in  the  Kato  and  Phillips  tank. 

A suggestion  was  made  by  Turner  (1973)  that  the  entrainment  velocity  should 
be  expressed  in  terms  of  and  a velocity  and  length,  ux , , characteristic  of 

the  rms  velocity  and  the  integral  length  scale  in  the  mixed  layer.  Turner  suggested 
that  Uj.  ~fS  and  D and  (1)  implies 

u« 


3 

-e  ~ R’  2 
Ul 


(2) 


where  Rix  is  the  turbulent  Richardson  number  expressed  in  terms  of  th,  ux 
and  j lL.  Long  (1975),  however,  presented  a uniform  theory  for  all  cases  in  which 
the  buoyancy  flux  near  the  interfacial  layer  q0  ~u8/D  where  u2  is  the  order  of  the 
velocity  in  the  mixed  layer  near  the  interface.  This  implies  an  Rig1  law  for  u« 
if  u2  is  the  velocity  used  to  scale  u,  and  in  the  expression  for  Ri8.  Long  suggested 
that  the  two  results  for  the  experiment  without  shear  could  be  reconciled  by  assuming 
that  the  small  density  variation  in  the  mixed  layer  was  sufficient  to  change  the  order 
of  magnitude  of  the  turbulent  velocities  from  ux  ~ fS  in  the  upper  part  of  the  layer 

1,  4_ 

to  Ug  ~ fS  Ri  6 or  Ug“P.  This  was  investigated  experimentally  by  Hopfinger 
and  Toly  who  concluded  that  the  velocities  were  everywhere  proportional  to  f and 
that  the  Rix  2 law  is  correct.  It  is  important  for  later  purposes  to  notice  that  their 
rms  velocity  measurements  were  not  made  particularly  close  to  the  interface  so 
that  the  observations  cannot  be  said  to  disprove  the  theory  of  the  author  that 


I 


q0  ~ u|/D.  To  avoid  confusion  between  the  paper  by  Long  (1975)  and  the  present 

1 

paper,  we  remark  that  the  present  paper  again  finds  q0  ~u|  /D,  but  that  ua  * f . 


r i 
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2.  Governing  Parameters. 

We  begin  by  considering  the  best  choice  of  parameters  governing  the  phenomena 
in  the  experiment.  The  external  quantities  are  the  frequency  of  oscillation  of  the 
grid  f,  the  stroke  S,  lengths  Mj, , Ma , . . . describing;  the  geometry  and  position  of 
the  grid,  and  the  initial  density  variation.  We  assume  that  the  dimensions  of 
the  tank  are  large  enough  to  be  neglected.  In  a recent  paper  (Long,  1977)  the 
author  has  shown  that  the  grid  may  be  replaced  by  a virtual  source  of  energy  at 
a horizontal  plane.  The  action  of  the  source  is  determined  by  a single  parameter 
K having  the  dimensions  of  viscosity  and  proportional  to  the  constant  eddy  viscosity 
in  the  turbulent  fluid  above  the  source.  When  stratification  exists,  the  eddy  viscosity 
will  be  constant  sufficiently  close  to  the  energy  source  since  the  velocities  are 
very  high  there  and  buoyancy  effects  negligible  as  we  will  see.  The  integral  length 
scale  is  proportional  to  the  distance  z from  the  virtual  source  so  that  if  ux  is  of 
the  order  of  u in  the  lower  portion  of  the  mixed  layer  at  some  height  zx , we  may 
replace  K by  u^j,  when  discussing  the  relation  of  experiments  to  theory. 

The  role  of  the  density  stratification  may  be  examined  by  integrating  the 
equation 


9q  _ 86 
3z  ~ 9t 


(3) 


first  over  the  mixed  layer  and  then  over  the  IL.  Let  us  assume  a linear  buoyancy 
field  in  the  non-turbulent  upper  layer  with  buoyancy  gradient  N2.  The  mean 
buoyancy  in  the  mixed  layer  is  nearly  constant  with  height  and  equal  to 

b„  = b - N2  (D+  h)  + £b  (4) 


where  b is  constant.  We  get 


=DT  -NaD-3I<D+h) 


(5) 
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If  the  buoyancy  field  at  the  initial  instant  was  the  linear  field  b^-  N2  z,  Eq.  (8) 
becomes 

> h'  N2  a 

EH- | Zb  =|  (D+h)2  (9) 

and  N is  the  fundamental  parameter  characterizing  the  buoyancy  field.  Otherwise 
the  constant 

V2  = ip+l/b  - ^(D+h)2  (10) 

and  N may  be  used  as  fundamental  constants.  If  the  upper  fluid  is  homogeneous, 
V3  = D+  j / Si  is  constant.  Using  h = aD,  we  obtain  DZb  = v2  is  constant  and 
v is  the  basic  parameter  characterizing  the  buoyancy  effect. 


3.  The  Interfacial  Layer. 

The  energy  equation  at  any  level  z is 

0 = - f(2p/p0+ua+va+w^  -wb  -e  (11) 


where  the  first  term  is  the  energy  flux  divergence,  u,  v,  w are  the  velocities, 
p is  pressure,  p0  is  the  reference  density,  q = -wb  is  the  buoyancy  flux  and  s 
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is  the  energy  dissipation.  We  must  have  turbulence,  although  perhaps  intermittent, 
in  the  IL.  Since  the  only  buoyancy  flux  occurs  in  the  turbulent  regions,  we  get  at 
any  level 

q = -Aj,  Ugbal  (12) 

where  Ax  is  a universal  constant,  u3  is  the  order  of  the  velocity,  b3  is  the  rms 
buoyancy  fluctuation  and  I is  the  intermittency  factor.  The  turbulence  is  certainly 
strongly  influenced  by  buoyancy  in  this  layer  and  it  is  reasonable  to  assume  that 
kinetic  and  available  potential  energies  are  of  the  same  order  not  only  in  the  waves 
but  in  the  turbulent  patches,  i.  e.  , 

(13) 

where  6 is  of  the  order  of  the  dimensions  of  the  patch,  and  u3  is  the  rins  velocity 
in  the  IL.  We  get 

(14) 


2 ,,  CD 

U3  ~ 6b3  ~ 6 -r- 


6 A Uo 
U = 


(hZ!b) 


Using  (13),  equation  (12)  becomes 


a _ _ A jgym 

q ‘ Ajh  u. 


(15) 


where  A2  and  A3  are  universal  constants.  Let  us  now  compute  the  dissipation. 
This  occurs  only  in  the  turbulent  patches  so  that 


ul  _ u| . (h&l 

6 h u3 


(16) 


so  that  c ~q.  Since  they  are  both  dissipative,  it  follows  that  they  are  of  the 
order  of  the  energy  flux  divergence.  At  the  upper  boundary  of  the  IL,  the 
kinetic  energy  of  the  waves  has  been  so  reduced  by  losses  to  potential  energy 
and  dissipation,  that  there  can  no  longer  be  wave  breaking  and  turbulence. 


rr 
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Thus  h is  the  depth  of  penetration  of  the  turbulence.  When  the  upper  fluid  has 
a uniform  density,  the  motion  above  is  irrotational  and,  neglecting  the  very  small 
time  changes  of  kinetic  energy,  the  energy  flux  must  be  zero  just  above  z = D+h 
and  therefore  zero  just  below.  Then  the  increment  of  energy  flux  across  the  IL 
is  equal  to  the  energy  flux  into  the  bottom  of  the  layer.  If  there  is  a density 
gradient  in  the  fluid  above  the  IL,  energy  flux  may  exist  there  because  of  the 
existence  of  internal  gravity  waves.  At  the  height  z = D + h,  however,  the 
energy  flux  is  too  weak  to  support  turbulence  so  that  it  has  decreased  to  a 
value  well  below  that  at  the  bottom  of  the  IL.  Therefore,  in  all  cases,  the 
increment  in  energy  flux  is  proportional  to  the  value  at  the  bottom  of  the  IL. 

The  constant  of  proportionality  will  vary  with  the  density  variation  in  the  upper 
layer.  Integrating  between  levels  just  outside  of  the  layer,  we  get  for  the  mean 
buoyancy  flux 


qh~uf  (17) 

where  u2  is  the  rms  velocity  in  region  R2  just  outside  the  IL.  Obviously  since 
pressures  are  continuous,  there  can  be  no  change  of  order  of  magnitude1  of  the 
velocity  across  the  level  z = D so  that  ua  ~ us  at  z = D.  Using  Eq.  (7)  we  may 
write 

Ci  Ua  ^ dAb  , h dAb  , Ab  dh  , 1 4V_  dD 
* = D"dT+  3_dr+^dt+2^dt 

. rt  (18) 

1 More  carefully,  we  recognize  that  the  bottom  of  the  IL  is  in  wave  motion 
with  occasional  small  areas  of  breaking.  The  front  between  the  turbulence  below 
and  the  mostly  laminar  flow  in  the  IL  is  very  sharp.  Instantaneously  the  pressures 
on  this  surface  are  of  order  u|  on  the  lower  side  and  therefore  u2  on  the  upper  side. 
The  latter  corresponds  to  velocities  of  order  u2  ~ u3  in  the  lower  portion  of  the  IL. 
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H 2 d 
= -D^+  ND-^(h+D) 


(19) 

(20) 


where  I0  and  60  are  evaluated  at  a level  just  above  z = D and  cx  is  a constant 
of  order  one  which  may  vary  with  the  stability  of  the  upper  layer. 

It  is  of  interest  to  compute  the  energy  flux  in  the  lower  part  of  the  IL.  In 
the  eddies  the  contribution  is  of  order  Iu2.  Thus  the  transfer  is  due  primarily 
to  the  waves  and  is  of  order  of  u2p  ~ u2c  where  c is  the  speed  of  the  energy  con- 
taining  waves  of  length  X , i.  e. , c ~X(t-  J • Since  u2c  is  of  the  order  of  u| , 
c ~ us  or 


This  suggests  that  the  eddies  in  the  region  R2  of  length  X are  resonating  with 
these  waves.  Since  X is  of  the  order  of  both  wave  length  and  amplitude  of  the 
breaking  waves,  then  60  should  be  of  order  X.  Comparison  of  (19)  and  (21)  shows 
that  this  is  indeed  the  case.  The  bottom  surface  of  the  interfacial  layer  will  have 
waves  of  length  and  amplitude  60  so  that  the  layer  R2  may  be  considered  to  be 


! 


of  thickness  60 . 


4.  Turbulence  in  the  Mixed  Layer  and  Final  Results. 


So  far  we  have  simply  assumed  turbulence  below  the  IL  without  considering 
its  properties  in  detail  and  the  relation  of  these  to  the  flux  q0  or  to  the  entrainment 
velocity  u, , Let  us  consider  the  effect  of  buoyancy  in  the  R2  layer  where  the 
magnitude  of  u is  u2.  Here  q ~ u2b2  ~ . Thus  b2h  ~uf  so  that  the  ratio  of 


kinetic  to  potential  energy  of  the  eddies  is 

Up  h 

bs6  ~ 6 


This  is  very  large,  so  that  the  buoyancy  variation  is  unimportant  dynamically 
in  R2.  This  will  also  be  true  in  the  mixed  layer  as  a whole  so  that  in  the  whole 
mixed  layer  we  may  assume  that  the  energy  dissipation  function  ‘.s  determined 
by  factors  independent  of  the  buoyancy.  The  dissipation  function  at  any  height  z 
in  the  mixed  layer  must  depend  only  on  K,  z and  D,  i.  e. , u3  / 1 - f(K,  z,  D) 
where  i is  the  eddy  length.  Thus  u3  z4  / 4K3  = f(z/D).  At  height  z = D-60  ^ D, 
the  length  l is  of  order  60  and  the  velocity  is  us  so  that 

— 3^..  - A 


where  Au  is  a universal  constant.  Using  (18) -(20)  we  may  write 


h4D6(Ab)4 


cK*  _ ^ dAb  , h dAb  . Ab  dh  . i dD  . h.  d ...... 

TZ 4 = D“ST  + Q-S~+7r7ir+3  ^ (D+T)^-(D+h)  (25) 


3.  3 

io  _ . K2  JL  ,4 

D _A4  D3\Ab/ 


-AsIo  ^ = D ^ - N*D^-  (h+  D)  (27) 

where  A4  and  As  are  universal  constants  and  c is  a constant  of  order  one, 
varying  with  the  buoyancy  gradient  in  the  upper  layer. 


We  now  make  the  assumption  based  on  the  observations  discussed  in 
Section  1 that  the  thickness  of  the  interfacial  layer  is  proportional  to  D and 


is  independent  of  the  stability.  When  the  upper  fluid  is  homogeneous, DAb  = v 
is  constant  and  in  terms  of  the  fundamental  Richardson  number  Ri  = vaD2/K2 
we  obtain 


U„D  p T 

pT-  = C4R1 

(28) 

- -f 

Io  = 02  Ri  4 

(29) 

^0  - „ Rf  * 
D - CV3 K1 

(30) 

u5D 

(31) 

where  are  constants.  As  shown  by  the  author  (Long,  1977), ux  (or  K)  is 
proportional  to  the  grid  frequency  f so  that  Eq.  (28)  leads  to 


j^RT-*-  , Ri’  = ff.  (32) 

_3 

compared  to  Ri*  2 suggested  by  experimenters.  Inspection  of  their  data, 
however,  reveals  little  reason  to  choose  one  law  in  preference  to  another. 

Notice  that  Eq.  (28)  may  be  integrated  to  yield 

2 

D«  t1"37 

The  results  in  (28) -(31)  are  the  same  for  the  erosion  of  a linear  density  gradient 

where,  however,  the  constants  may  vary  with  the  stability  of  the  upper  layer. 

i 

The  entrainment  law  now  yields  D « t . 

5.  Summary  and  Conclusions. 

It  is  now  possible  to  form  a reasonable  description  of  the  state  of  affairs 


i*~r  p,*  _ DAb 


in  the  experiment.  The  oscillating  grid  is  nearly  equivalent  to  a source  of 
energy  on  a plane  (z  = 0)  as  discussed  by  the  author  (1977).  A single  parameter  K 
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characterizes  the  source  where  K has  dimensions  L2T_1.  If  the  fluid  is  homo- 
geneous, the  source  started  at  t = 0,  causes  a turbulent  layer  to  develop  of  depth 
D ~ (Kt)*  separated  by  a front  from  the  non-turbulent  fluid  above.  Thus  the 

_i 

front  propagates  at  a speed  proportional  to  t 2.  After  the  front  has  moved  far 
away,  conditions  ultimately  reach  a steady  state  with  energy  supplied  by  the 
source.  In  any  layer  there  is  a balance  between  the  energy  entering  at  the  lower 
plane  less  the  energy  leaving  through  the  upper  plane  and  the  energy  dissipation. 
The  rms  velocity  components  are  proportional  to  each  other  with  universal  con- 
stants of  proportionality.  The  horizontal  component,  for  example,  is  given  by 


nz  _ t>  K 
*K-Bl  v/ 


(33) 


and  the  integral  length  scale  by 


- = B8 . - ' 


l 

z 


(34) 


where  v is  the  molecular  viscosity.  The  energy  dissipation  e decreases  rapidly 
with  distance  from  the  grid  as  kVz4.  If  a passive  additive  is  present,  its  mean 
concentration  S is  determined  by  solutions  of  the  equation 


as 


8s  S 


at  Kh  8z 


(35) 


where  Kh  is  the  constant  eddy  diffusivity  proportional  to  K. 

If  the  fluid  is  initially  stratified,  the  energy  source  creates  a mixed  layer 
of  thickness  D,  increasing  much  more  slowly  with  time,  separated  from  the 
non-turbulent  fluid  above  by  an  interfacial  layer  of  thickness  h.  Although  h is 
small  compared  with  D,  it  remains  proportional  to  D.  The  mixed 
layer  has  a very  weak  buoyancy  gradient  so  that  dissipation  nearly  balances  the 
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energy  flux  divergence,  and  the  loss  of  kinetic  energy  to  potential  energy  is 
small.  This  is  reflected  in  the  fact  that  the  kinetic  energy  u2  is  much  greater 
than  the  available  potential  energy  Xb  where  X is  the  eddy  scale  and  b is  the 
perturbation  buoyancy.  The  eddy  length  X is  proportional  to  distance  from 
the  source,  becoming  of  order  D in  the  center  of  the  layer.  X decreases  again 

_3 

toward  the  interface  becoming  of  order  60  ~ DRi  4. 

The  rms  velocity  u varies  as  K/z  in  the  upper  part  of  the  mixed  layer, 
decreasing  with  z but  maintaining  a proportionality  to  K.  In  the  experiment,  if 
viscosity  is  negligible  K is  proportional  to  the  frequency  of  the  grid  f.  As  we 
approach  the  interface  the  eddy  size  decreases  ultimately  to  60  in  the  vicinity 
of  the  interface  which  is  disturbed  by  waves  of  length  and  amplitude  of  order  60 . 
Since  the  energy  dissipation  is  of  order  u3/X  , the  dissipation  tends  to  increase 
by  the  reduction  of  X but  this  is  offset  by  the  loss  of  kinetic  energy  to  heat  which 
reduces  u.  If  the  stability  is  infinitely  great,  i.  e.  (K2  /v^D*)  - 0 , the  rms 
velocity  tends  to  zero  at  the  interface.  For  finite  but  small  values  of  K/vD, 

„ _1  2 4 

we  have  u2  ~ (K/D)Ri  4.  This  indicates  u2  “f2  instead  of  fT  as  originally 
suggested  by  the  author.  The  difference  lies  in  the  mechanism  for  the  decrease 
of  u as  we  approach  the  interface.  The  author's  original  suggestion  was  that 
the  buoyancy  gradient,  although  weak  in  the  mixed  layer,  was  sufficiently  strong 
to  decrease  the  kinetic  energy  at  the  expense  of  potential  energy.  In  the  present 
theory  u decreases  because  of  dissipation  to  heat  as  the  eddy  size  decreases 
near  the  interface.  However,  u does  not  go  to  zero  but  to  a value  corresponding 
to  a small  but  non-zero  eddy  length  50,  proportional  to  Ri  *.  This  dimension 
is  also  the  dimension  of  the  layer  over  which  the  interface  moves  up  and  down. 
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This  length  was  measured  by  Hopfinger  and  Toly  as  proportional  to  Ri*-1 

A 

where  Ri*  is  proportional  to  Ri. 

The  intar  facial  layer  is  very  stable  and  has  intermittent  turbulence  with 
- -x 

intermittency  factor  I ~ Ri  . The  turbulence  in  patches  of  dimensions  60  is 
caused  by  breaking  waves  of  length  and  amplitude  6.  These  are  energized  by 
resonance  with  the  pressure  fluctuations  on  the  lower  interfacial  surface  of 

i 

frequency  us/60  equal  to  the  natural  frequency  N = (£b/h)2.  Equating  these 

_3 

and  putting  h ~ D leads  to  ^ ~ DRi  4 as  in  Eq.  (30). 

Finally  the  buoyancy  flux  varies  in  both  the  mixed  layer  and  the  interfacial 

layer.  It  is  a maximum  q0  at  the  level  z = D where  it  is  of  order  Ug/D  determined 

by  the  rms  velocity  in  the  mixed  layer  close  to  the  interface  as  suggested  by  the 

author  in  an  earlier  paper.  q0  ~ u„  Zt>  where  u,  is  the  entrainment  velocity  or 

rate  of  increase  of  D and  is  proportional  to  Ri  T . This  is  close  to  measurements, 

j. 

but  differs  slightly  from  the  Ri*  2 law  proposed  by  the  experimenters.  The  en- 

2 1 . 

trainment  law  leads  to  a variation  of  D as  t1"1  compared  to  D a P when  the  upper 
fluid  has  a linear  density  gradient.  In  the  case  of  a linear  gradient  in  the  upper 
fluid  the  results  are  as  stated  above  except  that  the  constants  of  proportionality 

i 

may  vary  with  the  stability  of  the  upper  fluid.  The  t^  law  reveals  a slower  rate 
of  increase  for  the  linear  case  as  expected  but  it  is  somewhat  slower  than  that 
measured  by  Linden  (1975). 
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Figure  1.  The  container  has  a fluid  with  variable  mean  buoyancy  profile.  The 
lower  layer  of  depth  D is  fully  turbulent  and  has  a weak  buoyancy  gradient. 
The  buoyancy  decreases  strongly  in  an  interfacial  layer  (IL)  of  thickness  h 
between  the  mixed  layer  and  the  non-turbulent  layer  above,  the  buoyancy 
difference  across  h being  Ab.  The  upper  layer  has  a buoyancy  gradient  b(z). 
The  turbulence  is  caused  by  a grid  oscillating  up  and  down  with  stroke  S and 
frequency  f.  The  geometry  and  location  of  the  grid  is  given  by  lengths  , 
M2 An  important  non-dimensional  number  is  Ri*  = DAb/f2  S2. 

Figure  2.  Kato  and  Phillips  experiment  with  mean  buoyancy  profile.  The  turbu- 
lence in  the  upper  layer  is  caused  by  a screen  moving  along  the  surface 
exerting  a stress  t = u2  , where  u*  is  the  friction  velocity.  An  important 
non-dimensional  number  is  Ri*  =DAb/u2  . 
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1-  SUPPLEMENT  ARY  NOTES 


12.  SPONSORING  MILI  TAR  Y ACTIVITY 

Office  of  Naval  Research 


A theory  is  developed  for  turbulence  in  a stably  stratified  fluid’,  for  example  in  the 
experiments  of  Rouse  and  Dodu  and  of  Turner.  In  these  there  is  no  shear  and  the  turbulence 
is  induced  by  a source  of  energy  near  the  lower  boundary  of  the  fluid.  <A  growing  mixed  layer 
of  thickness  D appears  in  the  lower  portion  of  the  fluid,  separated  from  the  non-turbulent  fluid 
above,  in  which  the  mean  buoyancy  gradient  is  given,  by  an  interfacial  layer  (IL)  of  thickness  h 
The  lower  mixed  layer  ha3  a very  weak  buoyancy  gradient  and  the  large  buoyancy  difference 
across  IL  is  tb.  As  indicated  by  the  experiments  of  Thompson  and  Turner  and  Hopfinger  and 
Toly,  and  derived  by  the  author  in  a recent  paper,  if  u is  the  rms  horizontal  velocity  and  l is 
the  integral  length  scale,  ujt  = K is  a constant  in  a layer  of  homogeneous  fluid  agitated  by  a 
grid  and  at  some  distance  from  the  grid.  ^>When  there  is  stratification,  the  fluid  motion  is 
unaffected  by  buoyancy  forces  in  the  mixed  layer,  so  that  ui.  should  also  be  constant  in  the  lower 
portions  of  the  mixed  layer.  Since  A is  proportional  to  distance,  we  may  conveniently  suppose 
that  the  source  of  the  disturbances  is  at  a level  zC=J)  where  u is  infinite  in  accordance  with 
uz  = K.  Thus  we  may  take  K to  be  a fundamental  parameter  characterizing  the  turbulent  energy 
source.  Then  z is  distance  above  the  plane  of  the  virtual  energy  source. ' If  the  upper  fluid  is 
of  uniform  buoyancy,  DAb  = v may  be  shown  to  be  constant  if  we  accept  the  experimental 
observation  that  h is  proportional  to  D.„  In  general  v may  be  taken'to  be  a fundamental  paramet 
expressing  the  stability.  The  quantity  Ri  =vsDa/Ka  is  the  most  fundamental  of  the  several 
Richardson  numbers  that  have  been  introduced  in  this  problem  because,  with  its  use,  '^constant 
of  proportionality  do  not  depend  on  the  molecular  coefficients  of  viscosity  or  diffusion  (for 
high  Reynolds  number  turbulence)  or  on  the  geometry  of  the  grid. 

/ J„.  c!  1 7 P/fc  5 ' < continued  over) 


FORM 

I NOV  «S 


1473 


(PAGE  I) 


tcurity  Classification 


Security  Classification 


> 4 

KEY  WORDS 

LINK  A 

LINK  B 

WT 

ROLE 

WT 

RO  l 

Turbulence  in  stable  fluids' 
Entrainment  of  stable  interface 
Entrainment  velocitiy 
Intermittency 

Resonance  with  pressure  fluctuations 

13.  Abstract  (continuation) 

a _ 7 

The  theory  contains  a number  of  results:  (1)  u,  D/K  ~ Ri  T , where  u,  is  entrainn 
velocity,  i.  e.  the  rate  of  increase  of  D.  This  implies  D <*  t^1  for  a homogeneous  upper 
and  D “ t^  for  the  upper  linear  density  field.  The  entrainment  law  compares  with  a -■§  1 
suggested  by  several  experimenters;  (2)  The  turbulent  field  in  IL  is  intermittent  with  int 

A A _3 

mittency  factor  I ~ Ri  . The  turbulent  patches  have  dimension  6 ~ DRi  * ; (3)  If  the  rir 
horizontal  velocity  at  some  fixed  level  not  far  from  the  grid  is  denoted  by  ult  we  find  thi 
rms  velocity  near  the  interface  is  u3  ~ Uj,  Ri  4 ; (4)  The  buoyancy  flux  near  the  interface 
may  be  expressed  as  qo  ~u|/D  as  suggested  by  the  author  in  an  earlier  paper. 

A recent  experimental  paper  by  Hopfinger  and  Toly  proposes  that  uk  ~ ua.  The  pr 
theory  suggests  that  the  measurements  of  the  rms  velocity  were  made  too  far  from  the  ii 
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